Astronomy & Astrophysics manuscript no. lblucy 


© ESO 2013 


January 1, 2013 





Bayesian inference for orbital eccentricities 

L.B.Lucy 



Astrophysics Group, Blackett Laboratory, Imperial College London, Prince Consort Road, London SW7 2AZ, UK 



(N 

o 

(N 
O 

Q 

m 
P? 

6 

o3 



> 

(N 

m 
o 

(N 



13 



Received ; Accepted 



ABSTRACT 



Highest posterior density intervals (HPDI's) are derived for the true eccentricities e of spectroscopic binaries with 
measured values e w 0. These yield upper limits when e is below the detection threshold eth and seamlessly transform 
to upper and lower bounds when e > e t h- In the main text, HPDI's are computed with an informative eccentricity 
prior representing orbital decay due to tidal dissipation. In an appendix, the corresponding HPDI's are computed with 
a uniform prior and are the basis for a revised version of the Lucy-Sweeney test, with the previous outcome e = now 
replaced by an upper limit exj. Sampling experiments with known prior confirm the validity of the HPDI's. 

Key words. Stars: binaries: spectroscopic - Methods: statistical - Methods: data analysis 



1. Introduction 

For over two centuries, astronomers have been able to de- 
tect and analyse orbital motions for objects beyond the 
solar system. From measured positions (visual binaries) or 
radial velocities (spectroscopic binaries, exoplanets), orbital 
elements and their standard errors are typically obtained by 
least-squares. Accordingly, at this late date, when orbital 
motion is detected, prior information concerning similar ob- 
jects is available and can be incorporated into the analysis. 

An example is the Lucy-Sweeney (1971; LS) test for 
the statistical significance of small measured eccentricities 
e for spectroscopic binaries (SB's). The prior information 
that provided support for the LS test was as follows: 

1) The small e's (typically ^ 0.05) of numerous catalogued 
SB's were ~ E(e\0), the expected value of e due to mea- 
surement errors when the true value is e = 0. 

2) Savedoff's (1951) investigation of ecosuj, where to is the 
longitude of periastron. If an SB is also an eclipsing binary 
(EB), e cos uj can be determined from the velocity curve 
(sp) and independently from the light curve (ph). For EB's 
with secondary eclipses midway between consecutive pri- 
mary eclipses, (e cos uj) p h = 0, but the values of (e cos ui) sp 
are scattered over the interval (-0.05, 0.05), confirming that 
non-zero e's & 0.05 are often spurious. 

3) Tidal dissipation gives e — > as t — > oo since, for two 
point masses, e = is the state of minimum orbital energy 
for fixed orbital angular momentum. Thus, if the time con- 
stant for this decay is short enough, the system is likely to 
be observed when e <C ^(e|0) - i.e., well below the mea- 
surement threshold. 

In view of this prior information, LS adopted e = as 
the preferred (null) hypothesis (Ho) and imposed a mod- 
erately demanding level of significance before rejecting Hq 
and accepting an elliptical orbit. 

Given that secular evolution due to tidal dissipation was 
already well established in 1971 and is not less so now, there 
is merit in explicitly incorporating this mechanism into the 



analysis rather than implicitly via the LS preference for 
e = 0. This can be achieved by replacing the jrequentist 
approach of LS by one based on Bayes' theorem. 

2. Estimating the true eccentricity e 

We suppose that radial velocities of a single-lined specto- 
scopic binary (SB1) have been analysed to estimate the 
orbital elements and their standard errors. We ask: what 
can be inferred about the error-free elements when orbital 
evolution due to tidal dissipation is taken into account? 



2.1. Posterior probability 

We adopt the notation used in Lucy(1974; L74). The vec- 
tors of the estimated and the error-free elements are de- 
noted by x and £, respectively; and the distribution of 
probability in a;-space for given £ is denoted by Tl(x\£)dx. 
Integrating over £-space, we find that the distribution of 
probability in cc-space is <p(x)dx, where 



(x) = J ^(€)n(x|€)d€ 



(i) 
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Here tp is the probability density function (pdf) that rep- 
resents our prior knowledge about ip(£)d£, the distribution 
of probability in £-space for the SBl's true elements. 

Since SBl's with small e's are of interest, we suppose 
that Sterne's (1941) elements have been chosen, as in LS. 
The six elements are then: P, the orbital period; 7, the 
systemic velocity; K, the semi-amplitude of the velocity 
curve; To, an epoch at which the mean longitude is zero; 
and the pair ecostu, esinu. In terms of these parameters, 
the radial velocity curve of an SB1 is given by Eq.(l) in LS. 

With this choice and the assumption e 2 -C 1, the 
off-diagonal elements of the least-squares matrix have 
zero expectation values when observational weight is uni- 
formly distributed in phase (LS). Since observers strive 
to meet this condition, we assume it to be true; and this 
then implies negligible correlations between the elements. 
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Accordingly, to a good approximation, the error-broadening 
kernel II(a:|£) is simply the product of the six gaussians giv- 
ing the independent error distributions of the six elements. 

Next consider the pdf ?/>(£), which is convolved in Eq.(l) 
with the kernel n(£e|£). It follows that ip is only rele- 
vant if it varies significantly within the error bars of an 
individual orbital element. This is not true for P, 7, K, 
and To. However, when orbit circularisation is taken into 
account, ip(£) may vary significantly within the domain 
(e cos u> ± a ecosuJ ,esinuj± cr es i nuJ ) . Accordingly, we now in- 
tegrate Eq.(l) with respect to the (P, 7, K, T )-components 
of x and use the normalization of U(x\^)dx with respect 
to each of these four elements. The result is 



4> c (x 



s) = / 



ipc(€, v)H-c{x -£,y- n)d^dn 



(2) 



where x = e cos oj and y = e sin oj are the two surviving 
estimated quantities, whose error-free values are £ = ecosw 
and n — e sinzu, respectively, and the subscript c indicates 
that the pdf's are defined on a Cartesian grid. 

We now assume that the standard errors of e cos u> and 
e sin oj are equal, which is true of their expectation values 
when observational weight is uniformly distributed in phase 
(LS). The broadening kernel is then the circular normal 
distribution 



n c (z-£,y-r?) 



1 



-exp 



2/z 2 



where [i = o~ ecosLJ = o~ esmu! . 

Bayes' theorem can now be invoked to derive - see 
Eq.(lO) in L74 - the posterior probability 

Thus, if the least-square solution is 

x = e cos uj ± u, y = esinoj±u (5) 

the distribution of probability in (e cos m,£ sin tzj)-space 
is given by Q c {^,n\x,y)d^dn. To evaluate this posterior 
distribution, we must specify what our expectations were 
for ip c (^,r])d^dr], the prior distribution of probability in 
(e cos m, e sin tz7)-space. 

For the problem at hand, polar (p) coordinates (e, m) 
are more convenient than the Cartesian coordinates (£, rf) 
in the above formulae. The Jacobian of the transformation 
is J = e, so that, by conservation of probability, 



and 



ip p (e,w) =eip c (£,ri) 
Q p (e, m\x, y)=e Q c (£, n\x, y) 



(6) 
(7) 



2.2. A physical model for ip p 

Given P, there is a 2-D family of binaries that match the 
measured K at some inclination. For simplicity, a repre- 
sentative example is chosen to avoid integrating over all 
possibilities. 

Now consider an ensemble of such binaries that form at 
a uniform rate in the solar neighbourhood and have lifetime 
t*. We further suppose that all have e = £q at t — and 



that thereafter e decays exponentially with fixed e-folding 
time t*/v, so that 



In e{t) = In e — v 



(8) 



It follows that lne is uniformly distributed in the interval 
(lne*, lne ), where £* = eoexp(— v). The probability that 
£ G (e, e + de) is therefore dejve. 

If we now assume randomly oriented orbits, the prob- 
ability that w G (vj,vj + dm) is dm/2n. Accordingly, the 
prior probability that a binary is in the element dedw at 
(e,vj) is tppdedtu, where 



il)p{e,m) 



2ttve 



for £ G (£*,£ ) 



(9) 



Note that ip p decreases with increasing v = ln(eo/e*). 
Nevertheless, normalization of this pdf is maintained by 
the corresponding decrease in £*, the lower limit for inte- 
grations over e. 

In introducing this physical model, we in effect adopt 
an informative prior. The following quote is apt: " The real 
power of Bayesian inference lies in its ability to incorporate 
'informative' prior information, not 'ignorance' " (Feldman 
& Cousins 1998). 



(3) 2.3. Distribition of £ 



For numerical calculations, it is convenient to transform 
the integral in Eq.(2) into an integration with respect to 
the polar coordinates £ and zu, so that 



. JO 



tp p (£, u7)U c (x — £, y — n) d£dw (10) 



and to use logarithmic spacing in £ in order to accurately 
evaluate the contribution near £* . Note that 4> c is indepen- 
dent of uj when tp p is independent of w, as in Eq.(9). 

With (j) c evaluated, the distribution of probability in 
(e, tn)-space is given by Q p (£, w\x, y)d£dvj, with Q p = £Q C 
from Eq.(4). This 2-D pdf, which in general is not indepen- 
dent of w, may be of interest when analysing a particular 
SB1. But here our interest is in £, so we integrate over vj 
to obtain 



q(e\e) 



2tt 



Q p (e,zu\x,y) dm 



(11) 



Note that q(£\e) is independent of oj because of the absence 
of a correlation term in II C - see Eq.(3). This in turn fol- 
lows from the assumptions (Sect. 2.1) that e 2 <C 1 and that 
observational weight is uniformly distributed in phase. 

The posterior probability that the true eccentricity G 
(e, £ + efe) is therefore q(£\e)d£, with mean value 



< £ > 



=1: 



£ q(£\e)d£ 



(12) 



2.4. Bayesian terminology 

In the above, notation and terminology is from L74. To 
modern Bayesians, n c (x — £,y — n) is the likelihood and 
4> c {x,y) is the Bayes' factor. Elsewhere, modern usage is 
followed with respect to the terms prior pdf, posterior pdf 
and credible intervals. 
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3. Numerical results 

The theory of Sect. 2 is now illustrated by computing a par- 
ticular case in detail. 

3.1. Parameters 

There are two basic parameters, e//x and v, the number of 
c-folding decay times in i*. 

Note that eo is not a consequential parameter provided 
that e <C eo- In effect, we assume that an SB1 with e ~ 
has reached this configuration due to secular evolution and 
not due to the formation mechanism. In these calculations, 
e = 0.5. 

We choose v — 8.52, so that £* = 10~ 4 . Then, with 
[i = 0.01, s(t) < 2A5fi, the LS threshold, when t > 0.35t*. 
Thus, from our ensemble of SBl's, « 65% would be assigned 
e = by the LS test. 

3.2. The posterior pdf x(l°g£|e) 

Because of the concentration of probability towards £*, 
plots are more informative if the abscissa is loge rather 
than e. Accordingly, we define 

X(loge|e) = eg(e|e)xlnlO (13) 

In Fig.l, this pdf is plotted for e/n = 1.0,2.45,3.03 and 
3.72, values selected as follows: If e = 0, LS showed that 
the probability p e of exceeding e is given by 

,„„« = -! (£f (14) 

provided that Therefore, when testing H , the above 

values of e/fi correspond to levels of significance 61, 5, 1 and 
0.1 %, respectively. The criterion e//i > 1 was proposed and 
implemented by Luyten (1936); the 5% level by LS. 

For e/fi — 3.72, \ i s an asymmetric bell-shaped func- 
tion peaking at ps e, but with a tail extending down to 
loge* = —4.0. As el [i decreases, the peak weakens and the 
tail strengthens. At e//i = 1.0, the peak is absent and all 
the probability is in the tail, which derives from the phys- 
ical model. Intermediate calculations show that the peak 
first appears at ej[i = 1.42. Thus, for e//i < 1.42, x is a 
monotonically decreasing function of e. For e//i > 1.42, x 
is unimodal. 

The pdf x for e/fi = 3.72 and fj, = 0.01 in Fig.l is 
computed with e^j \i — 50. Repeating this calculation shows 
that x is independent of the upper limit e provided that 
e /n ^ 10. 

3.3. Percentiles 

The posterior pdf's in Fig.l imply asymmetric and rapidly 
changing credible (or Bayesian confidence) intervals ase//j 
varies. These are plotted in Fig. 2 for the indicated values 
of the probability that the true eccentricity is < e. For 
the normal distribution, the values 0.159,0.500,0.841 and 
0.977 correspond to displacements of — 1,0, +1 and +2cr, 
respectively. 

Fig. 2 reveals a dramatic switch in solution regime at 
ej//i w 3.6. For e ^ e^, the measured value e is close to the 
50% percentile and is tightly enclosed by the ±1ct intervals. 
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Fig. 1. The pdf x(loge|e) for the values of e indicated by the 
vertical lines. These are at e/fi = 1.00, 2.45, 3.03 and 3.72, cor- 
responding to levels of significance of 61, 5, l, and 0.1%, respec- 
tively. The bold curve is the pdf at the point where the LS-test 
switches from accepting to rejecting a circular orbit. 
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Fig. 2. Credible intervals for the posterior pdf x(log e|e) as func- 
tions of e. The plotted boundaries correspond to the true eccen- 
tricity being < e with the indicated probabilities. The values of 
e corresponding to levels of significance 61% (Luyten) and 5% 
(LS) are indicated, as is the regime transition at e t . 



In this regime, inferences are dominated by the actual mea- 
surement e ± fi. But this ceases to be so for e <, e^. Thus, 
for e/fi < 2.65, the 'solution' e — e falls outside the ±lcr 
intervals. Evidently, for e ^ e-f, inferences are increasingly 
dominated by the model of Sect. 2. 2. 
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Fig. 3. Highest posterior density intervals for the posterior pdf 
x(loge|e). The enclosed posterior probability of each HPDI is 
95%, and the prior pdf ip is given by Eq.(7). The filled circles 
are the posterior means < e > computed from Eq.(10). 



3.4. Highest posterior density intervals 

For e//j, < 1.42, x is a monotonically decreasing function of 
e (Sect. 3.2). It follows that traditional, equal-tail credible 
intervals exclude the point (e = £*) with greatest probabil- 
ity density (pd). This undesirable feature is avoided by in- 
stead computing highest posterior density intervals (HPDI; 
Box & Tiao 1973). These intervals are such that every point 
included has a higher pd than every point excluded. 

In general, the calculation of HPDI's is non-trivial. But 
here the pdf's \ are n °t pathological (Fig. 1), and so the 
following clipping algorithm finds the HPDI for specified e 
and designated enclosed probability (1 — a): 

Let log £fc L , . . . , log Eku be consecutive grid points 
that belong to and define the HPDI {ek L ,£k v )- Then an 
HPDI with smaller included probability is obtained by 
eliminating the grid point log Sk L if Xfc L < Xk v or the 
grid point log eh v if Xk v < Xk L ■ This is repeated until the 
included probability = (1 — a). 

The 95% HPDI's thus obtained for x(loge|e) are plot- 
ted for e = 0.000(0.001)0.080 in Fig.3. For e < 0.018, 
the HPDI's are effectively one-tail intervals since the lower 
bound is sl = £* = 10~ 4 . Accordingly, for this problem, 
HPDI's provide a seamless transition from upper limits for 
non-detected to two-sided intervals for detected eccentric- 
ities (cf. Fcldman & Cousins 1998). This is an appealing 
aspect of HDPI's for interpreting measured e's and their 
uncertainties. 

For e < 0.0142, the interval excluded from an HPDI 
is a single-tail because of the aforementioned monotonic- 
ity. For e > 0.0142, the pdf x is unimodal (Fig.l), with 
a maximum whose location — > e as e/fj, increases. This 
emerging, measurement-driven maximum eventually brings 
about the transition from one- to two-tailed intervals. For 
e > 0.018, the excluded probability is contained in two 



tails, with the upper tail's probability being initially 5%, 
but this decreases to w 2.5% when e//U » 1 because of x' s 
increasing symmetry - see Fig.l. 

Note that for the normal distribution, Af(0, 1), the 95% 
HPDI is the familiar equal tail interval (—1.96, +1.96) and 
the width of this interval is the narrowest that ecloses 95% 
of the probability. For non-symmetric pdf's, HPDI's are the 
narrowest intervals enclosing probability (1 — a) and as such 
are a natural generalization of the conventional equal-tail 
intervals for symmetric, bell-shaped pdf's. 

Confidence intervals are an economical means of con- 
veying the compactness or otherwise of a variate's distri- 
bution. The resulting loss of information, if of concern, can 
be avoided by plotting the pdf's, as in Figs.l and A. 2. 

3.5. Detection threshold 

The transition from one- to two-tailed HPDI's is a natural 
definition of the detection threshold e t h for non-zero ec- 
centricity - i.e., the measured value e above which attribu- 
tion purely to measurement errors is implausible. However, 
Fig.3 shows that £l remains <C e for a considerable inter- 
val beyond eth = 0.018, which in any case depends on is, a 
parameter likely to be only crudely estimated. 

If detection is crucial for a subsequent investigation - 
e.g., an observing program - then a threshold closer to e^ 
should be adopted (Figs. 2 & 3). 

3.6. Simulation 

The role that HPDI's can play in reporting eccentricities is 
best illustrated by sampling experiments. Accordingly, syn- 
thetic data for the model of Sect. 2. 2 are created as follows: 
If z\, Z2 are random numbers in (0, 1), a random ensemble 
member in (e, nj)-space is at 

£ = £o exp{— vz{) m = 2ir z 2 (15) 

Then, if Ci> C2 are random gaussian variates, this ensemble 
member is observed at the point 

e cos lj — e cos w + e sin cu — e sin vj + (16) 

Repeated N times, the resulting e's comprise a simulated 
observing campaign of N random ensemble members whose 
exact eccentricities e are known. 

With £» = 10~ 4 , a sample of N = 1000 SBl's are plot- 
ted in Fig. 4. As expected, the large majority of the points 
fall within the 95% HPDI's (e L ,eu). With Luyten's crite- 
rion e > fi, 79.5% of this sample would have their elliptical 
orbits accepted. But Fig. 4 clearly shows that most systems 
with e G (0.01, 0.03) have £ < 0.01 and so exceed Luyten's 
criterion because of the bias of the non-negative estima- 
tor e (LS). With the LS criterion e > 2.45/i, the accepted 
percentage drops to 40.5, and most of the systems with 
e G (0.01,0.03) would now be assigned e — 0, a marked 
improvement. 

The further improvement provided by the HPDI's is 
that the assignment e = can now be replaced by an 
upper limit. Thus, for example, with this choice of prior, 
an SB1 with e = 0.01 ± 0.01 is preferably reported as 
£ < £95 = 0.014. From the standpoint of testing theories of 
tidal dissipation, an upper limit is more informative than 
£ = 0. 
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Log c 

Fig. 4. Ensemble simulation comprising N = 1000 SBl's with 
the parameters of Sect. 3.1. The plotted points randomly sample 
the time interval (0,t„), have true eccentricities e from Eq.(6), 
and have gaussian measurement errors added to derive e. The 
95% HPDI's (sl,su) are indicated. For e < e th = 0.018, the 
lower limit el = £* = 10~ 4 . 



The sampling procedure can also be used to validate the 
upper limits. A sample with N = 10 6 has 161,521 systems 
with loge 6 (—2.1, —1.9) and 8,317 of these have e > £95. 
Thus 94.85% lie below the upper limit, closely agreeing with 
the designated 95%. 

3. 7. Upper limits 

Upper limits ejj when e < eth are the most useful products 
of this Bayesian machinery. With their validity confirmed 
above, their dependence on the orbital decay rate is now 
explored. 

Consider the representative measurement e = 0.01 ± 
0.01, so that e/p = 1.0, well below the detection threshold. 
In Table 1, the HPDI 95% upper limits £95 are given for this 
e/p as a function of £*. The decrease of £95 with increasing 
v reflects statistical reality: if the e-folding time f*/i/Cf„ 
most systems will have e -C e and a correspondingly small 
£95. Note that £95 even drops below e when v <; 13. 

In analysing an observed system without a good esti- 
mate of v or, equivalently, of £*, the conservative approach 
is to suppose that £* is no more than a factor ~ 10 below 
the measured e, thus avoiding claiming too low an upper 
limit E\j. An even more conservative approach is to adopt 
a uniform prior - see Appendix. 

4. Conclusion 

In this paper, by incorporating a model of an SBl's secular 
evolution, Bayes' theorem is used to infer bounds on its ex- 
act eccentricity e given its measured value e±p. Because the 
system's lifetime t* is finite, the asymptote e — is never 
reached. Thus, in contrast to Luyten(1936) and LS, the 
statistical problem is not one of model selection. Systems 



Table 1. HPDI 95% upper limits £95 when e/p = 1.0. 



log £* 


u 


£95 / /i 


-3.0 


6.21 


1.71 


-4.0 


8.52 


1.37 


-5.0 


10.82 


1.15 


-6.0 


13.12 


0.98 


-7.0 


15.42 


0.85 


-8.0 


17.73 


0.74 


-9.0 


20.03 


0.65 



assigned e — by these earlier tests should preferably have 
upper limits ejj computed. 

As Fig. 4 shows, the Bayesian upper limits contain the 
systems for which e is significantly larger than e due to mea- 
surement errors and bias. Thus, a major historical cause of 
spurious e's is eliminated. But physical causes remain, such 
as those due to proximity effects or to additional line ab- 
sorption by gas streams. An example is £ TrA, for which 
Skuljan et al. (2004) improved the precision of the radial ve- 
locities by a remarkable factor of 100 and reported a small 
but highly significant e = 0.0140 ± 0.0002. However, the 
significant non-detection of the Keplerian third harmonic 
(Lucy 2005) invalidated this claim 1 . As precision improves, 
similar testing for the third and higher Keplerian harmonics 
is essential for confirming that an orbit is truly eccentric. In 
addition, an update of Savedoff's (1951) work would pro- 
vide numerous examples of spurious e's for investigation 
into physical causes other than measurement bias. 



Appendix A: Uniform prior 

The model of Sect. 2. 2 is not appropriate if the SB or star- 
planet system has additional components causing signifi- 
cant gravitational perturbations. In this circumstance, a 
sensible option is to assume a uniform prior for £, as is 
already common practice for exoplanets (e.g., Ford 2006; 
Eastman et al. 2012). Together with the assumption of ran- 
domly oriented orbits, the prior probability of the system 
being in dedw is then tppdedw , where 

tp p (e,zu) = ^- (A.l) 

which now replaces Eq.(9) in Sects. 2.3 & 3. 

From the resulting posterior pdf q(e\e), the 95% HPDI's 
and means < e > are plotted in Fig. A.l for e G (0.00, 0.08). 
Comparison with Fig. 3 shows that the HPDI's are nearly 
identical for e/p <; 5. However, for e/p ^ 2 - i.e., in the 
non-detection domain - the upper limits zjj in Fig. 3 are 
markedly lower, reflecting the effect of tidal circularisation 
in creating systems with e -C e. 

In contrast with Fig. 3, Fig. A.l shows that a uni- 
form prior results in a sharply defined detection threshold 
s-th/ 'fJ> = 2.49, which is gratifyingly close to the (frequentist) 
LS value e LS /p = 2.45 given by Eq.(14) for p e = 0.05. The 
thresholds for other critical levels are given in Table A.l, 
with the corresponding upper limits (bounds) z\j in Table 
A.2. 



1 Hearnshaw et al. (2012) have just reported eleven additional 
non-detections. 
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Fig. A.l. Highest posterior density intervals for the posterior 
pdf q(e\e). The enclosed posterior probability of each HPDI is 
95%, and the prior pdf ip is given by Eq.(A.l). The filled circles 
are the posterior means < e > computed from Eq.(10). 



40 



10 





/ \ CLtl 


//*= 


2.49 






fi = 


0.01 
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HFD1: £ < 0.0394 
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Fig. A. 2. Detection threshold. The posterior pdf q(e\e) for e t h = 
0.0249, the measured value for marginal detection (Sect. 3. 5). 
The hatched area with e > e v = 0.0394 contains 5% of the 
probability. Note that q(0\e t h) = q(su\eth)- 

For the representative measurement ej \x = 1.0 of Sect. 
3.7, the 95% upper limit from Table A.2 is 2.41/z. This ex- 
ceeds the corresponding values in Table 1, confirming that 
the uniform prior is the more conservative option. 

A.l. A revised Lucy-Sweeney test 

Even without evidence of additional components, an in- 
vestigator may be reluctant to base an analysis of orbital 
elements on uncertain estimates of tidal decay. If so, the 



Table A.l. Detection thresholds for e. 



a(%) 




eth/fJ- 


31.7 


1.52 


1.71 


10.0 


2.15 


2.21 


5.0 


2.45 


2.49 


1.0 


3.03 


3.06 


0.1 


3.72 


3.74 



Table A.2. Upper limits (bounds) for e. 



ejji £68.3/^ £9<)/i» £95/^ £99 £99.9/^ 



0.0 


1.00 


1.64 


1.96 


2.58 


3.29 


0.5 


1.06 


1.75 


2.08 


2.73 


3.48 


1.0 


1.28 


2.05 


2.41 


3.11 


3.89 


1.5 


1.66 


2.50 


2.88 


3.59 


4.38 


2.0 


(2.59) 


3.03 


3.41 


4.12 


4.90 


2.5 


(3.31) 


(3.87) 


(3.96) 


4.65 


5.43 


3.0 


(3.86) 


(4.52) 


(4.82) 


5.18 


5.96 


3.5 


(4.37) 


(5.04) 


(5.36) 


(5.97) 


6.48 


4.0 


(4.89) 


(5.54) 


(5.86) 


(6.49) 


(7.18) 



assumption of a uniform prior for e should be attractive. 
Physically, this corresponds to no secular evolution of e and 
a formation mechanism that uniformly populates the inter- 
val < e < 1 . Accordingly, a system with e w is assumed 
to have formed as such (cf. Sect. 3.1). 

This neutral standpoint is an attractive basis for a re- 
vised version of the LS test in which the previous accep- 
tance of a circular orbit (Hq) is now replaced by an upper 
limit. 

The revised LS test with 1 — a = 0.95 proceeds as 
follows: 

1) The eccentricity e ± fi is derived from the least 
squares solution. 

2) If e/fi > 2.49, this measured value e ± \i is accepted. 

3) However, if e//i < 2.49, the measured value is 
replaced by the upper limit £95 obtained by interpolation 
in column 4 of Table A.2. 

To illustrate this revised test, upper limits are given 
in Table A. 3 for six SBl's for which circular orbits are re- 
ported on the first page of Table 1 in LS. Thus for YZ Cas, 
the least squares value e = 0.004 was rejected by the LS test 
and £ = accepted. We now compute the 95% upper limit 
as follows: The estimate /i = 0.0037 is derived from Eq.(13) 
and Table 1 of LS. Linear interpolation in Table A.2 at 
e/fj, = 1.1 then gives £95/^ = 2.50, whence e < £95 = 0.009. 

Table A. 3. shows that upper limits can differ by large 
factors, reinforcing the earlier remark (Sect. 3.6) that upper 
limits are to be preferred in testing theories of tidal dissi- 
pation. A critical data base of detections and upper limits 
would facilitate progress in this field. 

Acknowledgements. I am grateful to the referee for pointing out an 
error in statistical terminology. 
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Lucy: Orbital eccentricities 



Table A. 3. Upper limits £95 for SB1 sample. 



Star 


e 


A* 


e/n 


£95/^ 


£ 


YZ Cas 


(0.004) 


0.0037 


1.10 


2.50 


< 0.009 


HD 7345 


(0.046) 


0.019 


2.47 


3.93 


< 0.073 


DM Per 


(0.09 ) 


0.074 


1.22 


2.62 


< 0.19 


HD 16589 


(0.008) 


0.0087 


0.92 


2.35 


< 0.021 


HD 18337 


(0.073) 


0.030 


2.48 


3.93 


< 0.12 


HD 21912 


(0.005) 


0.011 


0.45 


2.07 


< 0.023 
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